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Abstract 


This  study  investigated  the  use  of  the  program  PROFILE 
for  the  design  of  low-speed,  low-observable  configuration 
airfoils.  For  our  purposes,  low-observable  configuration  is 
defined  as  being  characterized  by  a  small  leading-edge 
radius  of  less  than  1.18  percent  chord  and  a  thickness-to- 
chord  ratio  of  2.5  to  5.5  percent.  A  methodology  was 
developed  whereby  the  input  parameters  to  prescribe  the 
velocity  distribution  over  the  airfoil  could  be  determined 
by  a  power  law  relationship.  This  relationship  enables  the 
designer  to  develop  symmetric  airfoils  with  the  desired 
thickness-to-chord  and  leading-edge  radius  within  the  stated 
constraints.  The  resulting  symmetric  airfoils  compared  very 
well  with  NACA  4-digit  airfoils  of  like  thickness-to-chord 
ratio.  The  resulting  airfoils  have  an  increased  stall  angle 
compared  to  the  NACA  4-digit  airfoils. 
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APPLICATION  OF  THE  PROGRAM  PROFILE  FOR  THE  DESIGN  OF  LOW- 


SPEED,  LOW-OBSERVABLE  CONFIGURATION  AIRFOILS 


I.  Introduction 


Background 

In  the  1950 's,  the  United  State's  reconnaissance  and 
intelligence  communities  made  the  first  serious  attempts  to 
make  less  detectable  military  aircraft.  Through  these 
efforts,  a  standard  for  detectability  was  established  to 
compare  and  measure  the  performance  of  any  system.  This 
standard  is  called  the  radar  cross-section,  (RCS),  which  is 
determined  by  first  measuring  or  calculating  the  amount  of 
radar  energy  reflected  from  the  target  to  an  observer. 

Then,  the  RCS  is  determined  by  the  size  of  a  reflective 
sphere  that  would  return  the  same  amount  of  radar  energy. 
From  this  point  in  time,  the  race  for  low-observability  has 
continued  and  has  expanded  from  the  reconnaissance  world  to 
all  military  aircraft. 

Most  recently,  the  Stealth  fighter's  success  in 
Operation  Desert  Storm  proved  the  military  advantage  and 
lethality  of  low-observcible  technology  for  today's  Air 
Force.  It  also  emphasized  the  importance  this  technology 
will  play  in  subsequent  tactical  aircraft. 

However,  a  problem  exists  for  today's  airfoil  designers 
dealing  with  this  technology.  Currently,  a  shortage  of  data 
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exists  in  the  open  literature  for  configurations 
representative  of  general,  low-observable  airfoils.  The 
problem  may  come  about  since  most  data  is  a  result  of  either 
reduced  radar  cross  section  airfoils  or  radar  absorbing 
materials  testing,  and  is  therefore  classified.  Beyond  the 
classification  problem,  the  data  is  usually  limited  to 
persons  with  special  access  clearances.  Therefore,  these 
reports  could  not  be  located  in  a  normal  classified 
literature  search. 

To  address  the  problem,  the  Flight  Dynamics  Directorate 
of  Wright  Laboratory  (WL/FIMM),  Wright-Patterson  Air  Force 
Base,  through  consultation  with  the  Naval  Air  Weapons 
Center,  Warminster,  Pennsylvania,  obtained  the  analysis 
program,  PROFILE.  This  program  was  developed  by  Richard 
Eppler  (1)  of  the  University  of  Stuttgart.  This  program 
applies  an  inverse  design  method  that  has  been  developed 
over  a  period  of  about  twenty  years.  An  airfoil  velocity  or 
pressure  distribution  is  prescribed  from  which  an  airfoil 
geometry  is  obtained.  The  program  consists  of  a  conformal 
mapping  method  for  the  design  of  airfoils  with  imposed 
velocity-distribution  characteristics,  a  panel  method  for 
the  analysis  of  the  potential  flow  about  a  given  airfoil, 
and  a  boundary-layer  method. 

Scope 

For  the  purpose  of  this  thesis,  we  will  define  the 
general  low-observable  airfoil  configuration  as  being 
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characterized  by  both  a  small,  leading-edge  radius  (r^)  and 
a  small  thickness-to-chord  ratio  (t/c).  More  specifically, 
for  this  investigation,  the  desired  airfoil  t/c  range  is 
2.5%  £  t/c  £  5.5%.  The  defining  factor  in  limiting  the 
leading-edge  radius  will  be  the  threat  frequency.  Better 
qualities  of  low-observability  are  obtained  if  the 
scattering  return  of  the  radar  energy  is  reduced.  A 
reduction  is  achieved  if  the  radius  of  the  illuminated 
surface  is  less  than  the  wavelength  of  the  threat  radar  (2). 
For  this  investigation,  a  threat  frequency  of  10  gigahertz 
was  specified  by  WL/FIMM.  This  results  in  a  wavelength  of 
three  centimeters  or  approximately  1.18  inches. 

The  investigation  is  based  on  the  assumption  that  the 
airfoil  has  a  minimum  chord,  c,  of  100  inches.  This 
specifies  an  upper  limit  of  0.0118  for  r^.  Note,  r^  has 
been  nondimensionalized  by  the  chord.  For  a  current 
perspective,  the  Advanced  Tactical  Fighter  has  a 
r^  «  0.0025. 

Objectives 

This  investigation  has  three  goals.  The  first  goal  is 
to  determine  the  applicability  and  define  any  limitations  of 
PROFILE  to  airfoils  with  small  leading-edge  radii  and  small 
thickness-to-chord  ratio.  The  second  goal  is  to  develop  a 
method  for  designing  airfoils  of  this  class  by  defining  a 
range  of  input  parameters.  The  final  goal  is  to  compare  the 
resulting  data  for  these  low-observable  airfoils  with 
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existing  NACA  4-digit  airfoil  data  of  equivalent  thickness- 
to-chord  ratio.  This  comparison  will  show  if  the  low- 
observable  airfoils  have  any  improved  aerodynamic 
characteristics  over  the  equivalent  thickness  NACA  4-digit 
airfoils. 
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Potential-Flow  Design  Method 

This  section  will  overview  the  basic  theory  used  in 
PROFILE.  Eppler  (1,3, 4, 6)  provides  a  more  detailed  analysis 
from  which  this  section  is  based.  The  design  method  is  based 
on  the  conformal  mapping  of  the  complex  C~plane  into  a 
complex  2-plane  (Figure  1).  In  the  C-plane,  there  is  an 
infinite  flow  around  a  unit  circle  given  by 

C  =  5  +  in.  (1) 

with  the  rear  stagnation  point  located  at  C  =  1/  the  complex 
potential  is  defined  by 

F(0  =  e--^«C  +  e-^°C'^  -  ~lnC  (2) 

2IT1 

and  the  circulation  given  by 

r  =  4rrCsin(a)  (3) 

where  C  =  1.  This  choice  of  C  provides  a  unit  flow  velocity 
at  infinity.  The  flow  is  mapped  into  the  z-plane  by  a 
complex  function  given  by 

(D 

2(C)  =  PiC  -  E  P-vC"'’  (Pi’‘0)  (4) 

v=0 


with  [ C I  ^1  for  a  convergent  power  series .  This  equation 
satisfies  the  following  conditions 


Z(oo) 


=  00 


(5) 


/  )  =1 

These  conditions  mandate  that  the  unit  circle  in  the  C”plane 
will  map  into  the  z-plane  as  an  infinite  flow  (C  =  1)  at  an 
angle  of  attack,  a,  around  an  arbitrary  airfoil  profile. 

From  complex  theory,  we  know  the  complex  conjugate  velocity 
vector  w  in  the  z-plane  is  given  by 


w  =  ve 


-id  ^  dF  ^ 
dz 


dF 

IL 

dz 

-dZ 


(6) 


where  v  is  the  modulus  and  0  is  its  argument.  From  Equation 
(2),  we  can  obtain  the  function  dF/dC.  Thus,  the  problem 
remaining  is  to  determine  z ( C)  so  that  along  the  airfoil 
profile  there  will  be  a  prescribed  velocity  distribution. 

However,  the  problem  is  defined  in  the  C-plane  on  the 
unit  circle  (  C  =  and  not  in  the  z-plane.  The 

variables  x  and  <p  are  related  by 


X  »  -|(1  +  cos  <p),  (7) 

where  c  is  the  airfoil  chord.  Since,  the  velocity  v(4>)  is 
specified,  the  problem  can  be  solved  by  taking  the  logarithm 
of  Equation  (6),  giving 

In-^  =  In^?  -  In^  =  -Inv  +  i0.  (8) 

dF  ^  aC 

This  equation  shows  that  -In  v  can  be  obtained  directly  as 
the  real  part  of  ln(dz/dF).  If  ln(dz/dF)  is  known,  then 
dz/dC  is  determined  by  Equations  (8)  and  (2),  and  the 
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mapping  function  z(C)  can  be  obtained.  Therefore,  a  trial 
solution  for  Jn(dz/dF)  is  introduced  in  accordance  with 
Equation  (8)  and,  with  several  manipulations,  an  equation  is 
introduced  as 

CD 

In^  =  In(l-r^)  ^  ^  (9) 

m=0 

The  function  dz/dC  resulting  from  Equation  (9)  must  be 
compared  to  the  requirements  from  Equation  (4)  and  must  also 
be  analytical  in  satisfying  the  condition  |C|  >1.  It  is 
important  to  note  that  Equation  (9)  can  be  ordered  in 
decreasing  powers  of  C •  A  comparison  of  this  form  of 
Equation  (9)  with  the  derivative  of  Equation  (4)  points  out 
that 

aQ  —  0  bQ  =  0  a^^  =  1  —  0.  (10) 

Examining  the  summation  in  Equation  (9)  along  the  boundary 
(C  =  e^*^),  the  real  and  imaginary  parts  can  be  expressed, 
respectively,  as  follows, 

CD 

P(<p)  =  5^  cos  mtp  +  sin  m<p)  (11) 

07=0 

CD 

Q((p)  =  53  (^0!  m<p  -  Bjj,  sin  m4>) .  (12) 

07=0 

Equation  (9)  can  thus  be  rewritten  as 
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dz 

‘dT 


-  ln(l  -  C'^) 


=  Pi<p)  +  iQ(<l>). 


(13) 


JC=e’^ 

Using  Equations  (2),  (6),  and  (9)  the  real  part  of  equation 
(13)  is  determined  to  be 


v(<p) 

P(<f>)  =  -In 

2 1  cos  ||  -  a  1 

(14) 


The  imaginary  part  is  the  conjugate  harmonic  function  of 
P((^),  which  is  defined  by 


2n 


QW  -  ^ 


P(i)r)  cot 


'\f  -  4> 


dijr. 


(15) 


In  this  method,  the  function  Q(<^)  is  numerically  evaluated. 
Now,  if  Equation  (14)  is  substituted  into  Equation  (13), 
after  some  manipulation,  it  is  determined  that 


(  — 1 

=  -4  sin  — 

1  cos  (—  -  I 

1  1 

\d<^/ 

2 

'  \2  )' 

v(0) 

Equation  (16)  can  be  numerically  integrated  and  then  split 
into  its  real  and  imaginary  parts  directly.  By  examination 
of  Equations  (13-16),  it  can  be  seen  that  the  airfoil 
coordinates  can  be  found  if  Q(^)  is  known. 

In  this  method,  the  velocity  distribution  is  not 
completely  arbitrary.  The  problem  is  still  constrained  by 
the  conditions  in  Equation  ( 5 )  and  those  subsequently  found 
in  Equation  (10).  With  this  information,  the  following 
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conditions  are  found  from  the  Fourier  expansion  of  Equation 

(11). 


27r 

I  P(<t>)d(f>  =  0 


77t 


^  P(<p)  cos  <p  d<p  =  1 


(17) 

(18) 


2n 

P(<p)  sin  <p  d(p  =  0.  (19) 

Therefore,  the  evaluation  of  the  airfoil  coordinates  is 
specified  by  Equations  (11-16)  under  the  conditions  of 
Equations  (17-19). 

In  this  method,  an  airfoil  is  produced  by  selecting 
angles  of  attack,  a  =  a*,  over  segments  of  the  unit  circle 
for  which  the  velocity  distribution  v{(p,a*)  is  to  be 
constant.  Thus  P{(p)  is  shown  to  be 

P(«)  .  -In _ _ 

2|cos  j|  -  a‘jl 

Equation  (9)  indicates  the  coefficients  Sjj,  and  bjj,  are  the 
keys  to  determining  the  mapping  function  z(C ) ,  which  gives 
the  airfoil  shape.  Thus,  Equation  (9)  indicates  the 
importance  of  P(<p)  in  the  determination  of  the  airfoil 
shape.  It  is  important  to  note  in  examining  Equation  (9), 
that  any  change  in  the  angle  of  attack,  a,  does  not  change 
the  mapping  function,  2(C).  Since  P(«^)  is  independent  of  a. 


(20) 
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the  velocity  distribution  for  any  arbitrary  a  is  found  from 
the  following  relationship: 


vl(p,a)  _  v*((p,a*) 


(21) 


The  invariance  with  respect  to  a  makes  it  possible  to 
designate  v*  by  choosing  values  of  a*  over  several  segments 
of  the  airfoil.  It  is  necessary  that  v(<t>,a)  be  continuous 
over  the  airfoil.  However,  this  does  not  constrain  a*  to  be 
continuous,  it  only  requires  that  v*((p,a)  counteract  any 
discontinuity  in  a*.  Thus,  the  following  matching  condition 
must  hold  between  the  values  of  a*  for  adjoining  arc 
segments: 


lim 

e=0 


(<Pi  * 
cos^_  - 


lim 

e=0 


(22) 


Examination  of  Equation  (20)  reveals  that  P{<f>)  is 
undefined  at  the  stagnation  point.  This  problem  is 
eliminated  by  introducing  an  arc  limit,  near  the 

leading-edge.  In  this  method,  is  not  specified  in  the 

input,  but  is  left  as  a  free  parameter  to  be  determined 
through  satisfaction  of  the  conditions  in  Equation  (10). 

The  unit  circle  is  divided  into  I  arcs,  where  each  arc 
i  extends  from  to  <p^,  with  the  following  definitions, 

(pQ  =  0  and  (p2  =  2n.  The  specified  velocity  distribution  is 
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represented  by 


v(<p,a*)  =  vj^  w(<t>)  ^  ‘P  ^  f  <23) 


where  is  a  constant  determined  from  solving  the 
transcendental  equation  for  ^i/xe  is  an  empirical 

function  introducing  the  main  pressure  recovery  and  a 
closure  contribution  on  each  surface  of  the  airfoil. 

A  discussion  on  Eppler's  (1,3)  meaning  of  main  pressure 
recovery  as  it  applies  to  this  c'ass  of  airfoils  will  be 
addressed  in  Chapter  III.  Eppler  gives  the  function  w(^)  as 


w{<p)  = 


\^^fcos  <t>  -  cos 

(cos  <b  -  cos  <ps]  2 

1  1  +  cos  <p„  j 

I  1  -  cos  03  j 

(24) 


for  (0  ^  <p  ^ 


w{<p)  = 


2 

1+K 

cos  0  -  cos 

► 

1-0.36' 

. 

cos  0  -  cos  03 

► 

1  +  cos  0^, 

1  -  cos  03 

(25) 


for  (0i,ie  <  ^  ^  271).  The  terms  enclosed  by  the  braces,  { 
},  are  to  be  dealt  with  as  a  specially  defined  function, 
f(<p) : 


0 


{f{<P))  = 


^  0 

f{<p)  >  0- 


(26) 


The  first  factor  in  both  Equations  (24)  and  (25) 
represents  the  main  pressure  recovery  contribution.  Figure 
3  shows  the  total  amount  of  pressure  recovery,  o,  which  is 
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defined  as  the  ratio  of  the  velocity  at  the  airfoil's 
trailing-edge  to  the  velocity  at  the  specified  recovery 
location.  The  length  of  the  recovery  is  specified  by  <f>^  and 
for  the  upper  and  lower  surfaces,  respectively.  We  see 
from  Equation  (24)  that  for  the  upper  surface  the  total 
cunount  of  pressure  recovery  is  represented  by 


^  v(0,a*)  ^ 

v{<p„,a*)  (l+cos<^^J 


(27) 


The  relative  slope  at  the  beginning  of  the  pressure  recovery 
((p  =  (p^)  is 


g/  = 


-  (1  +  cos  <p„) 


(28) 


As  we  can  see  from  Equations  (27)  and  (28),  the  total  amount 
of  pressure  recovery  and  the  shape  of  the  recovery  are  found 
by  specifying  any  two  of  w,  w'  ,  K,  or  jU. 

The  second  term  in  Equations  (24)  and  (25)  represent 
the  closure  contribution  for  the  trailing-edge.  Eppler  (4) 
suggests  this  area  is  usually  confined  to  a  small  area 
around  the  trailing-edge  (^g  «  24®  -  36°).  For  the  purposes 
of  this  investigation,  <pg  will  be  held  constant  at  24°  for 
all  the  subsequent  runs  for  both  the  upper  and  the  lower 
surface.  The  values  of  and  are  not  specified  in  the 
input.  They  are  determined  by  an  iteration  procedure  that 
ensures  that  the  profile  obtained  will  have  a  closed 
trailing-edge.  In  the  program,  Eppler  provides  nine 
trailing-edge  iteration  mode  options  (1,7).  This 
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investigation  addresses  symmetrical  airfoils  which  warrants 
choosing  iteration  mode  three  for  our  computations.  With 
iteration  mode  three,  all  upper-surface  are  replaced  by 
+  Aa  and  all  the  lower  surface  a^*  are  replaced  by  - 
Aa. 

In  summary  the  complete  specification  of  v(<p)  needs  the 
following  information: 

(1)  Arc  limits,  <p^  (not  including 

(2)  Angle  of  attack  values,  a^*  for  all  arcs 

(3)  Pressure  recovery  specifications 

(4)  Length  of  the  closure  contributions 

Thus,  we  are  left  with  I  constants  v^,  and  as 

mentioned  above,  and  Kjj  as  free  parameters.  Also,  there 
are  the  I  matching  conditions  from  Equation  (26)  and 
Equations  (17-19)  to  be  satisfied.  These  conditions  are 
satisfied  in  closed  form,  leading  to  a  transcendental 
equation  for  This  equation  is  solved  (Regula  falsi), 

allowing  for  the  determination  of  v^,  and  Kj^. 

Potential-Flow  Analysis  Method 

In  PROFILE,  Eppler  utilizes  a  panel  method  with 
vorticities  distributed  parabolically  along  each  panel  (1). 
With  the  airfoil  coordinates  determined  from  the  potential- 
flow  design  method,  the  geometry  of  each  panel  is  determined 
by  a  spline  fit  of  the  coordinates.  The  panels  are  defined 
with  the  airfoil  coordinates  as  their  endpoints.  At  each 
airfoil  coordinate,  the  condition  that  the  inner  tangential 
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velocity  be  zero  is  satisfied.  This  mathematical  condition 
corresponds  to  the  physical  condition  of  the  flow  being 
tangent  to  the  surface.  The  analysis  is  made  at  a  =  0“  and 
a  =  90°  and  any  other  angle  of  attack  is  derived  by  the 
superposition  of  these  two  solutions.  The  trailing-edge 
singularities  resulting  from  the  circulation  r  around  the 
airfoil  are  addressed  according  to  whether  the  shape  of  the 
trailing-edge  is  sharp  or  blunt  (1). 

The  importance  of  using  the  parabolic  vorticity 
distribution  for  this  application  is  shown  by  Hess  (5), 
since  this  class  of  airfoil  will  result  in  very  thin 
trailing-edges.  Hess  (5)  indicates  the  method  will  provide 
remarkably  accurate  solutions  for  airfoils  with  these 
characteristics . 

Boundary-Layer  Method 

Eppler  (1,4,6)  employs  an  integral  method  for  the 
analysis  of  the  boundary- layer .  A  basic  overview  of  the 
method  will  be  presented  here  for  clarity, 

Nondimensional  Variables.  In  this  method,  the  unit 
system  selected  is  in  terms  of  a  characteristic  length,  L,  a 
characteristic  velocity,  Vj^,  and  a  characteristic  density, 
pp,  of  the  problem.  In  this  system,  for  example,  time  is 
measured  in  L/V^  and  mass  in  PrL^.  As  a  result  of  this 
system,  in  an  incompressible  problem  the  only  density  in  the 
problem  is  the  reference  density,  or  p  h  i.  All  eguations 
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are  transformed  to  the  nondimesional  variables  by  setting 
Pj^  H  1,  Vp  =  1,  and  LSI.  For  the  problem  of  flow  around 
an  airfoil,  the  length  unit,  L,  is  the  chord,  c,  the 
velocity  unit,  Vj^,  is  the  free  stream  velocity,  U„,  and  the 
density  unit  is  Pr  =  Pto- 


General  Equations.  The  conventional  coordinate  system 
is  used,  where  x  is  streamwise  along  the  surface  and  y  is 
normal  to  the  surface.  The  governing  equations  for  this 
method  are  the  two-dimensional,  incompressible  boundary- 
layer  equations. 


u 


(1^)  * = -n  ^  I? 


(29) 


=  0 
-Sy 


(30) 


The  potential  flow  velocity  U(x)  is  defined  by 


U{x)  =  lim  u(x,y) . 

y-HD 


(31) 


The  displacement  thickness,  the  momentvun  thickness,  and  the 
energy  thickness  are  defined  respectively,  as  follows: 


03 


u(x,y) 

U(x) 


(32) 


62(x) 


CP 


r  u(x,y)l^  _ 

u(x,y)  \ 

i  U(x)  \ 

U(x)  } 

dy 


(33) 
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(34) 


63(X)  = 


0 

r  u(x,y) 

1  _  /u(x,y)\2 

\  U(X)  1 

dy. 


These  thicknesses  will  be  used  to  characterize  the  viscous 
effects  of  the  boundary-layer.  From  these  thicknesses,  the 
shape  factors  11^2  and  H32  are  defined  as 


Hi2 


(35) 


«32 


63 


(36) 


The  boundary- layer  equations  can  be  manipulated  and 
integrated  in  the  y  direction  with  the  integral  parameters 
above  to  obtain  an  integral  momentum  equation 

6'  .  («,  .  26,  )^  =  C,  .  ^  (37) 

and  an  integral  energy  equation 

6^  3  63  ^=C„*^.  (38) 

In  the  equations  above,  v^  is  the  normal  velocity  at  the 
surface,  where  blowing  is  defined  as  positive  and  suction  is 
negative,  and  the  primes  denote  differentiation  with  respect 
to  X.  No  blowing  will  be  considered  in  this  investigation, 
thus  Vq  =  0.  The  skin  friction  coefficient  Cf  is  defined  as 


Cf  = 


pU^ 


(39) 


and  the  dissipation  coefficient,  Cp,  is 
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2 

pU^ 


(40) 


Cd 


tD 


dy. 


Laminar  Boundary-Layers.  At  this  point,  the  method 
follows  Von  Karman's  idea  to  approximate  solutions  for 
laminar  boundary-layers.  It  is  assumed  that  u(x,y)  is  to 
be  approximated  by 


u(^,y)  (41) 

where  6(x)  and  A.(x)  are  unknowns.  The  value  of  X  determines 
the  velocity  profile  shape  with  respect  to  y/6,  and  6  is 
only  able  to  stretch  the  profile  in  the  y-direction.  Eppler 
(6)  discusses  how  Equation  (38),  used  under  the  assumption 
of  Equation  (41),  is  the  best  approach  since  the  shape 
factors  H22  and  H32  are  only  functions  of  X.  It  is 
important  to  note  that  H22  is  only  a  function  of  H32.  Now, 

Cf  and  Cj)  are  found  to  be 


C/ 


^*(^32) 

ReU62 


(42) 


2D*(H32) 

ReU62 


I7<oL  U62 

where  Ne*,  =  -  — t  =  ReU&j 

*2  V  2 


(43) 


based  on  the  local  flow  velocity  and  the  momentum  thickness 
in  accordance  with  the  unit  convention  described  earlier. 

As  is  apparent,  the  precision  of  the  approximate  method 
is  dependent  on  how  well  the  exact  solution  of  Equations 
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(29)  and  (30)  are  represented  by  Equation  (41).  Eppler  (6) 
indicates  that  the  Hartree  profiles  give  the  best 
approximation  in  regions  of  adverse  pressure  gradient 
(1.51509  <  H32  <  1.57258).  A  negative  dp/dx  is  defined  as  a 
favorable  pressure  gradient  and  a  positive  dp/dx  is  defined 
as  an  adverse  pressure  gradient.  The  value  of  H32  =  1.51509 
specifies  laminar  separation.  With  all  the  assumptions  made 
so  far,  the  functions  fli2i  £*/  and  D*  in  terms  of  H32  are 
derived  as 

if  H32  <  1.57258 

Hi,  =  4.02922  -(583.60182  -  724.5591^3,  +  227.18220^3,)* 

^  (44) 

U32  -  1.51509  ) 

e*  =  2.512589  -  I.686O96H32  +  0.39154h|2  -  0.031729h|2 

if  H32  >  1.57258 

Hi2  =  79.870845  -  89.58214^32  +  25.715786h|2  (45) 

e*  =  1.372391  -  4.226253^32  +  2.221687h|2 

D*  =  7.853976  -  10.260551^32  3.418898h|2- 

With  these  functions,  Equations  (37)  and  (38)  become  coupled 
ordinary  differential  equations,  which  can  be  numerically 
integrated. 

Turbulent  Boundary-Layers.  A  turbulent  boundary-layer 
can  be  addressed  with  Equations  (37)  and  (38),  also. 

However,  The  Reynolds  stresses  impart  new  problems  which  do 
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not  allow  the  functions  to  be  derived  as  efficiently  as  in 
the  laminar  case.  Eppler  (6)  provides  expressions  derived 
from  both  empirical  and  semi-empirical  investigations  of 
turbulent  boundary- layers  from  several  researchers.  The 
expressions  are 


_  IIH32  ^  15 
48^32  -  59 


(47) 


Cf  =  0.045716[(Hi2  “  (48) 

-1 

=  0.0100[(Hi2  ' 

With  these  expressions,  Equations  (37)  and  (38)  again  become 
coupled  ordinary  differential  equations  to  be  numerically 
integrated.  For  this  case,  it  is  not  possible  to  define  a 
fixed  value  of  H32  as  the  separation  point.  But,  it  is 
known  that  for  H32  >  1.58  there  will  be  no  separation  and 
for  H32  <  1.46  there  definitely  will  be  separation.  Eppler 
(4)  states  that  this  method  gives  lower  values  of  H32  for 
adverse  pressure  gradients  than  other  methods,  therefore 
turbulent  separation  is  assumed  to  occur  at  H32  =  1.46. 


Boundary-Layer  Transition.  The  transition  criterion 
for  this  program  is  based  upon  local,  boundary-layer  data 
(4).  The  expression  used  to  indicate  transition  is 
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(50) 

In  ^  18.4  H32  -  21.74  +  125^(^32  -  1.573)^  -  0.36r 

where  r  represents  a  roughness  factor  for  0  ^  r  i  6  and  m 
allows  switching  between  the  normal  criterion  (m=0)  and  a 
modified  criterion  (m=l).  For  the  roughness  factor,  r  =  0 
corresponds  to  smooth  conditions  and  r  =  6  corresponds  to  a 
very  rough  surface  or  a  very  turbulent  free  stream.  For 
this  investigation,  smooth  conditions  were  assumed.  The 
modified  criterion  was  developed  from  flight  test  data  to 
improve  the  performance  of  the  normal  criterion  for  cases 
where  the  shape  factor  was  near  laminar  separation,  at  high 
Reynold's  number  f  r  in  favorable  pressure  gradients.  For 
airfoils  in  ta'  -  class,  the  normal  criterion  and  modified 
criterion  results  were  not  significant,  so  all  data 
presented  is  with  the  normal  criterion.  When  transition  is 
indicated,  the  computations  switch  from  Equations  (44) 
through  (46)  to  Equations  (47)  through  (49).  Physically, 
transition  does  not  occur  at  one  point,  so  the  criterion  in 
Equation  (50)  indicates  where  the  laminar  boundary-layer 
ends. 
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III.  Analysis  and  Limitations 


For  an  inviscid,  incompressible  fluid,  Bernoulli's 
equation  can  be  written  as 

1  O 

p  +  '2^'^  ~  constant.  (51) 

It  follows  from  this  equation  that  in  a  given  flow  field,  an 
increase  in  velocity  results  in  a  reduction  of  pressure  and 
vice  versa.  In  incompressible  flow,  Bernoulli's  equation 
gives  the  relationship  between  the  local  values  of  p  and  V 
with  those  in  the  free  stream  as 


Po,  +  =  p  *  (52) 

which  can  be  rearranged  to 

P  -  Poo  =  (53) 

At  this  point,  it  is  useful  to  introduce  a  nondimensional 
difference  between  the  local  and  free  stream  pressure,  Cp. 
This  pressure  coefficient  is  defined  by 


Thus,  from  Equations  (53)  and  (54),  we  see  that 


(54) 


(55) 


The  local  V  is  composed  of  three  elements:  the  velocity 
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induced  by  the  thickness  of  a  symmetrical  airfoil  at  zero 
angle  of  attack,  the  velocity  corresponding  to  the  camber, 
and  the  velocity  corresponding  to  the  change  in  angle  of 
attack.  For  symmetrical  airfoils,  the  velocity  will  only  be 
dependent  on  the  thickness  and  change  in  angle  of  attack. 
Since  the  term  (V/V„)^  in  Equation  (55)  can  never  be 
negative,  it  follows  that  the  value  of  Cp  will  be  negative 
over  the  suction  side  of  the  airfoil.  This  will  occur  on 
the  upper  surface  of  the  airfoil,  since  the  flow  is 
accelerating  in  this  region.  Thus,  at  some  point  the  flow 
will  reach  a  minimum  pressure,  Cp^^„,  and  from  this  point  the 
pressure  recovery  will  start. 

Amount  and  Location  of  the  Pressure  Recovery 

Eppler  (4)  states  that  in  designing  airfoils  there  is 
no  general  rule  for  picking  the  parameters  of  the  recovery 
function.  Thus,  it  is  essential  to  define  the  parameters 
best  suited  to  develop  a  specific  class  of  airfoils. 

In  the  discussion  of  the  input  parameters,  only  the 
upper  surface  variables  will  be  specified.  All  equations 
are  identical  for  the  lower  surface,  with  the  only 
difference  being  that  the  lower  surface  parameters  are 
contrasted  from  the  upper  by  an  overbar. 

In  the  input  files,  <{>^  is  specified  by  the  parameter  k. 
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The  two  parameters  are  related  by  the  following  expression: 

^  (56) 

where  N  is  the  total  number  of  points  representing  the 
airfoil  minus  one  (the  point  at  the  trailing-edge  is 
included  for  both  the  upper  and  lower  surfaces).  Unless 
otherwise  noted,  the  value  of  N  =  60  will  be  used. 

In  the  previous  chapter,  the  function  w(<^),  which  was 
introduced  in  Equations  (25)  and  (26),  contains  a  factor 
that  Eppler  (1,4,6)  refers  to  as  the  main  pressure  recovery 
contribution.  Equations  (27)  and  (28)  defined  the  total 
amount  of  pressure  recovery,  w,  and  the  relative  slope  of 
the  recovery,  o'.  In  both  these  equations,  the  value  of  </>„, 
the  location  of  the  beginning  of  the  main  pressure  recovery, 
is  significant. 

First,  Eppler 's  terminology  for  the  pressure  recovery 
will  be  explored.  For  comparison  with  a  thin  airfoil,  the 
pressure  recovery  of  a  thick,  cambered  airfoil  will  be 
examined.  Figure  4  shows  pressure  distributions  for  Airfoil 
E398  (t/c  =  .1427),  a  thick,  cambered  airfoil  developed  by 
Eppler  (4).  For  this  airfoil,  the  main  pressure  recovery  is 
specified  to  begin  at  approximately  the  mid-chord  location. 
Figure  4  shows  that  through  a  =  6®,  where  a  is  relative  to 
the  chord  line  for  this  cambered  airfoil,  Eppler 's  term  of 
main  pressure  recovery  describes  very  well  the  location 
specified  by  <p^.  However,  for  a  >  6°,  a  pressure  recovery 
is  also  beginning  near  the  leading-edge  of  the  airfoil.  By 
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a  =  9°,  the  pressure  recovery  at  the  leading-edge  becomes 
significant  and  the  distinction  of  main  pressure  recovery 
becomes  blurred. 

For  the  class  of  thin  airfoils  being  examined,  Eppler's 
choice  of  terms  is  especially  confusing.  If  an  airfoil's 
t/c  is  small,  its  leading-edge  radius  will  also  be  small. 
Small  leading-edge  radii  on  low-speed  airfoils  may  cause 
early  separation,  or  leading-edge  stall,  as  the  angle  of 
attack  is  increased.  As  the  angle  of  attack  increases,  the 
magnitude  of  the  upper  surface  suction  peak  increases  and 
moves  forward.  Figure  5  shows,  for  Airfoil  G48 
(t/c  =  .0533),  that  the  dominant  pressure  recovery  occurs 
near  the  leading-edge  even  for  small  angles  of  attack  on 
thin  airfoils.  Figures  4  and  5  indicate  that  this  occurs  on 
both  thick  and  thin  airfoils,  but  the  effect  is  more 
pronounced  at  smaller  angles  of  attack  for  thin  airfoils. 
Thus,  we  see  it  is  not  necessarily  the  main  pressure 
recovery  that  is  being  defined  by  but  what  I  will  term  a 
secondary  pressure  recovery.  Figure  6  summarizes  the  terms 
used  in  the  description  of  the  pressure  distributions  used 
in  this  investigations. 

Figure  5  also  shows  the  influence  of  <p^  on  the  pressure 
distribution.  Moving  the  secondary  pressure  recovery 
(Eppler's  main  pressure  recovery)  from  the  trailing-edge  to 
the  leading-edge,  by  changing  the  value  of  X  in  the  input 
file,  results  in  a  damping  effect  on  the  dominant  pressure 
recovery  near  the  leading-edge.  The  important  question  to 
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answer,  in  obtaining  a  method  for  the  design  of  these 
airfoils  with  small  t/c,  is  defining  the  effect  of  the 
secondary  pressure  recovery  placement.  Figure  7  illustrates 
the  effect  of  the  placement  of  the  secondary  pressure 
location  on  t/c  and  It  is  apparent  that  the  choice  of 

X  =  0  or  29.99  (<f>^  =  0  or  »  180  degrees)  gives  the  smallest 
t/c . 

From  Equations  (27)  and  (28),  if  is  chosen  to  be 
zero,  by  definition  o  =  1  and  o'  =0.  Therefore,  K  =  0  and 
=  1  are  chosen  to  satisfy  both  equations.  Essentially, 
the  effect  of  this  choice  is  to  disregard  the  secondary 
pressure  recovery.  These  results  make  intuitive  sense, 
since  the  thickness  of  the  airfoil  is  related  to  the  amount 
of  pressure  recovery.  Figure  5  also  shows  a  higher  suction 
peak  near  the  leading-edge  for  X  =  0  or  29.99,  which  would 
indicate  an  airfoil  with  a  small  t/c  and  leading-edge 
radius.  Thus,  the  airfoils  designed  in  this  study  will  use 
a  value  of  =  A  =  0. 

Bouudary-Laver  Displacement  Iteration 

The  program  contains  the  option  to  add  the  displacement 
thickness  to  the  airfoil  shape  obtained  from  the  potential 
flow  design  (7).  This  will  allow  the  analysis  to  be 
performed  on  an  effective  airfoil  shape  in  a  potential  flow. 
This  displacement  iteration  is  performed  for  angles  of 
attack  specified  in  the  potential  analysis  section. 

The  displacement  thickness  is  computed  in  the  boundary- 
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layer  method.  The  displacement  thickness  is  added  to  the 
contour  for  the  specific  angles  of  attack  and  Reynolds 
numbers  specified.  The  resulting  shape  is  then  rerun 
through  the  panel  method. 

Unfortunately,  problems  were  encountered  with  the 
utilization  of  this  option  in  the  program  as  described 
above.  The  problem  arose  because  the  program  was  set  up  to 
only  output  the  new  data  through  a  plotting  subroutine. 

Also,  these  subroutines  were  written  for  the  plotting 
software  used  at  the  University  of  Stuttgart.  At  first, 
this  appeared  to  be  a  very  minor  inconvenience.  A  modified 
procedure  would  be  to  have  the  program  output  the  new 
velocity  or  pressure  distributions  before  the  data  was  sent 
to  the  plotting  subroutines.  However,  upon  close 
examination  of  this  section  of  the  code,  it  was  found  that 
the  pressure  or  velocity  distributions  for  the  new  shape 
were  recalculated  in  a  plotting  subroutine.  Even  with  this 
surprising  discovery,  it  appeared  the  information  could 
still  be  readily  retrieved.  Figure  8  shows  the  results. 
Unfortunately,  the  velocity  distribution  does  not  match 
Eppler's  data  (7).  After  many  attempts  and  much  time  spent 
to  find  the  error  in  the  approach,  it  was  concluded  that 
some  transformation  must  exist  within  the  plotting 
subroutines  that  was  not  accounted  for. 

Therefore,  an  input  file  was  created  with  the  new 
airfoil  coordinates,  which  were  obtained  in  the  main  section 
of  the  boundary- layer  calculations.  Figure  8  shows  that  the 
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new  airfoil  profile  matched  with  Eppler's  profile  (7).  The 
program  was  reexecuted  to  perform  the  potential  analysis  on 
the  airfoil  created  with  the  new  coordinates.  Figure  9 
shows  that  the  results  compare  very  well  to  those  of  Eppler 

(7) .  Figure  9  shows,  for  this  particular  airfoil  at  the 
specified  angle  of  attack  and  Reynolds  number,  there  is  not 
a  significant  difference  between  the  two  calculations. 

Further  confidence  was  gained  in  this  modified 
procedure  through  the  comparison  of  results  found  in  Graham 

(8) .  Figure  10  shows  a  significant  difference  between  the 
experimental  data  versus  the  potential  flow  solution  with  a 
compressibility  correction  for  M  =  .4.  But,  Figure  10  also 
shows  the  results  are  very  good  when  the  boundary-layer 
displacement  iteration  is  performed,  as  modified  above. 

The  determination  needed  to  be  made  to  see  if  this 
calculation  was  essential  to  the  analysis  for  the  thin 
airfoils.  The  boundary-layer  displacement  effects  on  are 
evaluated  in  a  very  simple  way  if  the  displacement  iteration 
is  not  used.  The  lift  curve  slope,  2ttt|,  is  reduced  to  2n, 
where  the  efficiency,  t^,  is  defined  by  Eppler  (4)  as, 

n  =  1  +  0.78t/c.  (57) 

This  reduction  in  lift-curve  slope  assumes  that  the 
potential  flow  thickness  effects  are  canceled  by  the 
boundary- layer  displacement  effects.  From  inviscid,  thin 
airfoil  theory,  where  the  potential  flow  streamline  is  at 
the  airfoil  surface,  the  lift-curve  slope  for  an  infinitely 
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thin  airfoil  is  2n.  For  an  airfoil  with  finite  thickness, 
Equation  (57)  corrects  for  the  resulting  increase  in  the 
lift-curve  slope  as  the  thickness  of  the  airfoil  increases. 
For  a  viscous  flow,  the  potential  flow  streamlines  are 
displaced  by  the  boundary-layer  thickness.  As  the  flow 
proceeds  along  the  airfoil,  the  boundary-layer  thickness 
increases.  This  results  in  a  reduction  in  the  lift-curve 
slope.  Their  combined  effect  proves  to  be  negligible, 
resulting  in  a  lift  curve  slope  very  close  to  the 
theoretical  value  of  2n  (9).  Figure  11  shows  that  the 
calculated  with  this  simple  method  compares  very  well  to 
experimental  data  for  a  NACA  0006  airfoil  at  Re  =  3.0  x  10^ 
(10).  Thus,  for  the  purpose  of  investigating  the  input 
parameters  necessary  to  define  a  method  of  obtaining  an 
airfoil  of  this  class,  it  was  not  essential  to  run  a 
boundary-layer  displacement  iteration. 

Viscous  Corrections  to  Section  Lift  Coefficient  (Separation) 
An  inviscid  section  lift  coefficient,  is  easily 
obtained  by 

Cl  =  2rrnsin(a)  (58) 

where  t)  is  the  efficiency  defined  by  Equation  (57).  Of 
course,  boundary-layer  separation  will  have  a  significant 
effect  on  the  Eppler  (4)  corrects  using  a  method 

derived  from  Helmholtz  potential  flow-wake  theory. 
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Considering  the  upper  surface,  the  correction  is  made  by 

ACj  =  2TrAa  (59) 

where 

In  Equation  (60),  s^^p  is  the  length  of  separation 
calculated  in  the  boundary-layer  analysis,  6yg  is  the 
airfoil  slope  near  the  trailing-edge,  and  is  the  angle  of 
attack  relative  to  the  chord  line.  Equations  (59)  and  (60) 
basically  result  in  a  reduction  of  the  angle  of  attack,  a, 
by  a  factor  Aa.  This  is  illustrated  in  Figure  12.  This 
same  procedure  is  used  for  the  lower  surface  with  the 
corresponding  parameters.  It  should  be  noted  that  the 
method  constrains  the  correction  to  be  negative  for  the 
upper  surface  and  positive  for  the  lower  surface. 

As  discussed  earlier,  an  airfoil  with  a  small  t/c  and 
small  leading-edge  radius  is  prone  to  early  separation.  As 
an  example.  Figures  13  and  14  show  that  Airfoil  G31  is  a 
symmetric  airfoil  with  t/c  =  0.0495  and  practically  no 
leading-edge  radius.  Figure  15  shows  Cj^(a)  for  this  airfoil 
at  Re  =  3.0  x  10®.  This  figure  shows  the  airfoil  stalling 
where  a  is  as  low  as  0.5®.  Here  we  see  a  problem  arising 
from  the  correction  defined  by  Equations  (59)  and  (60). 

After  the  stall,  the  program  is  predicting  a  negative  lift 
for  positive  a  between  approximately  0.5  and  3  degrees. 

This  is  not  a  realistic  result  for  a  symmetrical  airfoil. 
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Examining  Equation  (60),  it  is  seen  that  Sggp  is  the 
dominant  parameter  in  the  determination  of  Aa,  since  the 
trailing-edge  slope  for  thin  airfoils  are  on  the  order  of  a 
couple  of  degrees  or  less.  Hence,  this  viscous  correction 
will  predict  negative  values  of  after  the  occurrence  of 
stall  for  thin  airfoils  that  have  the  characteristic  of 
separation  moving  very  abruptly  toward  the  leading-edge  at 
small  angles  of  attack.  Looking  back  to  Figure  11,  we  see 
an  example  of  an  airfoil  where  stall  does  not  occur  at  a 
very  low  angle  of  attack.  The  NACA  0006  has  a  much  larger 
leading-edge  radius  than  Airfoil  G31  and  is  slightly 
thicker.  However,  the  stall  occurs  abruptly  for  these 
airfoils.  Again,  the  program  does  not  predict  correctly 
after  the  stall,  but,  since  the  angle  of  attack  is  greater, 
a  negative  is  not  predicted.  However,  it  is  important  to 
note  that  the  program's  results  compare  very  well  in  the 
calculation  of  the  magnitude  of  The  value  of 

will  be  used  in  the  comparison  with  NACA  4-digit  airfoils. 

Specification  of  the  Velocity  Distribution 

Presently,  two  of  the  four  requirements  for  the 
complete  specification  of  v(<p)  have  been  established,  the 
length  of  the  closure  contributions  and  the  secondary 
pressure  recovery  specifications.  Now,  the  determination  of 
the  arc  limits,  <p^,  and  the  values  for  the  corresponding 
arcs  must  be  made  such  that  it  results  in  an  airfoil  with 
2.5%  s  t/c  s  5.5%  and  a  small  leading-edge  radius. 
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As  pointed  out  in  Chapter  II,  the  function  P(0)  plays  a 
important  part  in  the  determination  of  the  airfoil  shape. 
Also,  the  invariant  nature  of  P(^),  with  respect  to  the 
angle  of  attack,  gives  us  the  relationship  between  the 
specified  and  the  resulting  velocity  distribution. 

This  can  be  seen  by  rearranging  Equation  (21),  resulting  in 


v{(f>,a) 


(61) 


Applying  Equation  (7),  Eppler  (3)  points  out  that  if  V(x,a) 
is  considered  with  a  certain  a  for  a  segment  on  the  upper 
surface  with  and  V(x,a*)  that 

^(V(x,a))>0  if  a  <  a*  (62) 

and 


(V{x,a)  )  <  0  if  a  > 


(63) 


Differentiating  Equation  (51)  with  respect  to  x  yields 


(64) 


which  gives  the  familiar  relationship  between  velocity  and 
pressure  gradient.  Comparing  Equations  (62)  and  (63)  with 
Equation  (64)  gives  us  a  very  important  insight  into  the 
effect  of  the  selection  of  ,  with  respect  to  the  a  of  the 
incoming  flow,  on  the  velocity  or  pressure  gradient 
produced.  When  a  <  a^* ,  a  favorable  pressure  gradient  will 
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occur  over  that  upper  surface  segment  and  if  a  >  a^* ,  an 
adverse  pressure  gradient  occurs. 

For  example,  in  the  previous  section  the  Airfoil  G31 
was  introduced  in  the  discussion  of  early  leading-edge 
separation.  Airfoil  G31  was  defined  by  two  arc  segments  on 
the  upper  surface  as  follows, 

=  I*  for  30  V  >  22.5 

=  0*  for  22.5  ^  V  ^  0  (65) 

where 

and  with  the  negative  values  of  a^*  for  the  corresponding 
segments  on  the  lower  surface.  Obviously,  these  input 
parameters  produced  an  airfoil  with  small  t/c  and  a 
practically  sharp  leading-edge.  However,  since  the 
magnitudes  of  the  are  small,  an  adverse  pressure 
gradient  develops  at  very  small  angles  of  attack,  which 
induces  a  rapid  separation.  Figure  15  shows  the  stall  for 
Airfoil  G31  occurring  at  approximately  a  =  0.5®. 

This  is  not  a  desirable  characteristic  for  low-speed 
airfoils.  What  happens  if  we  increase  the  value  of  on 

this  airfoil?  Figures  16  and  17  indicate  that  as  is 

increased,  both  t/c  and  leading-edge  radius  will  increase. 

On  the  upper  surface,  the  favorable  pressure  gradient,  or  an 
accelerating  velocity,  over  a  wider  a  regime  on  the  front 
portion  of  the  airfoil  results  in  a  thicker,  more  curved 
shape  at  the  leading-edge.  This  causes  a  paradox.  The 


32 


airfoil  needs  to  have  a  small  t/c  and  leading-edge  radius 
(small  values  of  a*),  yet  stalling  at  very  low  angles  of 
attack  is  undesirable. 

Figure  18  compares  the  pressure  distributions  of  a  NACA 
0006  and  a  NACA  0012  at  a  =  3®  and  a  =  8“  computed  with  the 
program.  Considering  the  upper  surface  or  the  suction  side 
of  the  airfoil,  Figure  18  illustrates  that  for  the  thinner 
airfoil,  the  favorable  pressure  gradient  region  is  very 
short  and  the  initial  pressure  recovery  has  a  steeper 
adverse  pressure  gradient  than  that  of  the  thicker  airfoil. 

The  previous  analysis  indicates  that  a  function  which 
produced  a  decreasing  a*  would  appear  to  be  a  good 
compromise  in  the  design  process.  Figure  19  shows  the  hand- 
drawn  curves  of  decreasing  a*  created  to  test  the  above 
analysis.  The  values  of  a*  from  the  TA2-TA5  curves 
(depicted  in  Figure  19  by  symbols)  were  input  into  the 
program  with  the  corresponding  arc  segments.  The  TA2-TA5 
curves  are  not  smooth  because  the  curves  shown  are  created 
with  only  the  fourteen  piecewise  continuous  values  of  a*, 
for  each  surface,  chosen  for  input  into  the  program.  The 
(Vi,aj^*)  pairs  for  the  curves  TA2-TA5  are  shown  in  Appendix 
A,  Table  A. 1. 

The  program  allows  for  as  many  as  28  arc  segments  to  be 
defined  for  one  airfoil,  including  both  upper  and  lower 
surfaces.  A  large  percentage  of  the  14  arc  segments  per 
surface  were  placed  near  the  leading-edge  of  the  airfoil  to 
capture  the  steep  slopes  of  the  functions  near  the  leading- 
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edge.  Figures  20-31  show  the  airfoil  profile,  leading-edge 
profile,  and  pressure  distributions  for  the  TA2-TA5  curves. 
Appendix  B  contains  an  introduction  to  the  significant  input 
lines  used  during  this  investigation.  Appendix  C  contains 
the  input  files  for  each  airfoil.  The  value  of  r^.  was 
estimated  from  the  figures  of  the  airfoil's  leading-edge. 
Table  1  summarizes  the  t/c  and  r^  for  the  TA-Series 
airfoils. 


TABLE  1 

SUMMARY  FOR  TA-SERIES  AIRFOILS 


Airfoil 

t/c  (%) 

■t 

Comments 

TA2 

12.05 

0.0170 

No  separation  thru  10° 

TA3 

9.05 

0.0158 

No  separation  thru  10° 

TA, 

5.56 

0.0100 

separation  between  8-10° 

TA5 

2.37 

0.0042 

separation  between  4-6° 

From  this  comparison,  the  TA4  and  TA5  curves  are  shown 
to  result  in  airfoils  near  the  desired  t/c  band  and 
prescribed  r^  limit.  These  curves  produce  airfoils  that 
fall  within  the  prescribed  r^  limit  and  are  just  outside  the 
desired  t/c  band.  However,  TA4  has  the  preferred  separation 
characteristic  compared  to  TA5. 

The  TA4  and  TA5  curves  shed  light  on  the 
characteristics  of  the  a*  distribution  needed  to  reach  the 
desired  objective,  since  these  curves  resulted  in  airfoils 
within  the  prescribed  t/c  and  r^  bands.  The  task  will  be  to 
determine  a  reliable  way  to  develop  the  input  parameters 
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needed  to  produce  an  airfoil  with  a  desired  t/c  and/or  rf 
To  achieve  this  goal,  the  best  approach  seemed  to  be 
that  of  choosing  a  mathematical  function  for  a*  which 
produced  the  characteristics  of  TA4  and  TA5  curves.  The 
above  analysis  indicates  that  curves  like  TA4  and  TA5  could 
be  represented  by  a  power  law  relationship  of  the  form 

a*  =  (x/c)"^  (0<d<l).  (66) 

This  relationship  will  give  a  large  value  of  a*  near  the 
leading-edge  that  will  decrease  rapidly.  Thus,  an 
investigation  of  the  applicability  of  this  function  for 
developing  this  class  of  airfoils  needed  to  be  made.  The 
program  will  not  accept  any  value  of  a*  ^  90°.  Imposing 
this  limit  on  Equation  (66)  constrains  the  value  of  the 
exponent,  d,  in  the  power  law  relationship,  such  that 
d  <,  0.592.  However,  it  is  obvious  that  values  of  d  >  0.592 
can  be  used  if  a  coefficient,  A  <  1,  is  included  in  the 
power  law  relationship  such  that, 

a*  =  A(x/c)~^  (0<d<l).  (67) 

Figures  32-36  show  the  a*  distribution  obtained  using 
Equation  (67)  for  the  following  five  values  of  A:  1,  1/2, 
1/3,  1/4,  and  1/5.  Equation  (67)  was  computed  with  a 
A(x/c)  =  0.0005.  This  spacing  is  needed  to  break  out  the 
arc  segments  required  to  define  the  leading-edge,  as  was 
done  for  the  TA-series  airfoils.  The  (v^jQ^^*)  pairs  derived 
from  these  power  law  curves  are  shown  in  Appendix  A,  Tables 
A. 2  -  A. 6 . 


35 


IV.  Results  and  Comparisons 


Case  Siunmary 

Figures  32-36  show  the  results  from  Equation  (67)  being 
evaluated  at  the  different  values  of  d  for  each  of  the  five 
cases  for  A.  Several  d  values  were  evaluated  for  each  case. 
The  values  of  d  were  determined  in  each  case  so  to  define  a 
clear  relationship  between  the  thickness-to-chord  ratio  and 
the  leading-edge  radius  throughout  the  specified  thickness- 
to-chord  regime.  The  establishment  of  this  relationship 
will  give  the  designer  the  ability  to  determine  the  input 
parameters  needed  to  achieve  a  certain  thickness-to-chord  or 
leading-edge  radius  requirement.  The  values  of  a*  are  taken 
from  the  curves  found  in  Figures  32-36,  in  the  same  manner 
as  done  for  the  TA-series.  The  values  for  the 
pairs  are  shown  in  Appendix  A,  Tables  A.2-A.6.  The  airfoils 
created  with  these  a*  distributions  will  be  referred  to  as 
the  T-series  airfoils.  In  the  following  cases  only  the 
airfoils  falling  into  the  specified  thickness-to-chord  and 
leading-edge  radius  bands  will  be  presented. 

Case  1.  A  =  1.  For  the  values  of  d  evaluated  at  this 
value  of  A,  four  airfoils  meet  the  specified  objectives  for 
the  thickness-to-chord  ratio  and  the  leading-edge  radius. 

The  airfoil  profile,  leading-edge,  and  section  lift 
coefficient  for  Airfoils  T25,  T30,  T3,  and  T4  are  shown  in 
Figures  37,  38,  40;  41,  42,  44;  45,  46,  48;  and  49,  50,  52, 
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respectively.  These  figures  show  that  as  the  value  of  d  is 
increased  both  the  thickness-to-chord  ratio,  the  I aadxug- 
edge  radius,  and  the  maximum  lift  coefficient  increase. 
These  airfoils  are  summarized  in  Table  2. 


TABLE  2 

T-SERIES  AIRFOIL  SUMMARY  (A  =  1) 


Airfoil 

d 

t/c  (%) 

ft 

Ci^ax  /  Stall  Angle 
(Re  =  3  X  10^) 

T25 

0.25 

2.96 

0.0020 

0.495  /  4.5® 

T30 

0.30 

3.59 

0.0030 

0.66  /  6.0° 

T3 

0.333 

4.09 

0.0040 

0.77  /  7.0® 

T4 

0.4 

5.17 

0.0075 

0.935  /  8.5® 

Case  2.  A  =  1/2.  For  the  values  of  d  evaluated  at  this 
value  of  A,  two  airfoils  meet  the  specified  objectives  for 
the  thickness-to-chord  ratio  and  the  leading-edge  radius. 

The  airfoil  profile,  leading-edge,  and  section  lift 
coefficient  for  Airfoils  T402  and  T22  are  shown  in  Figures 
53,  54,  56;  and  57,  58,  60,  respectively.  These  figures 
show  that  as  the  value  of  d  is  increased  both  the  thickness- 
to-chord  ratio,  the  leading-edge  radius,  and  the  meuEimum 
lift  coefficient  increase.  These  airfoils  are  summarized  in 
Table  3. 
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TABLE  3 


T-SERIES  AIRFOIL  SUMMARY  (A  =  1/2) 


Airfoil 

d 

t/c  (%) 

T402 

0.4 

3.44 

0.0038 

0.66  /  6.0® 

T22 

0.5 

4.75 

0.0060 

0.88  /  8.0® 

Case  3.  A  =  1/3.  For  the  values  of  d  evaluated  at  this 
value  of  A,  two  airfoils  meet  the  specified  objectives  for 
the  thickness-to-chord  ratio  and  the  leading-edge  radius. 

The  airfoil  profile,  leading-edge,  and  section  lift 
coefficient  for  Airfoils  T503  and  T603  are  shown  in  Figures 
61,  62,  64;  and  65,  66,  68,  respectively.  These  figures 
show  that  as  the  value  of  d  is  increased  both  the  thickness- 
to-chord  ratio,  the  leading-edge  radius,  and  the  maximum 
lift  coefficient  increase.  These  airfoils  are  summarized  in 
Table  4. 


TABLE  4 

T-SERIES  AIRFOIL  SUMMARY  (A  =  1/3) 


Airfoil 

d 

t/c  (%) 

ft 

/  Stall  Angle 
(Re  =  3  X  10®^) 

T503 

0.5 

3.79 

0.0050 

0.77  /  7.0® 

T603 

0.6 

5.08 

0.0085 

0.95  /  8.7® 

Case  4.  A  =  1/4.  For  the  values  of  d  evaluated  at  this 
value  of  A,  three  airfoils  meet  the  specified  objectives  for 
the  thickness-to-chord  ratio  and  the  leading-edge  radius. 

The  airfoil  profile,  leading-edge,  and  section  lift 
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coefficient  for  Airfoils  T504,  T604  and  T234  are  shown  in 
Figures  69,  70,  72;  73,  74,  76;  and  77,  78,  80, 
respectively.  These  figures  show  that  as  the  value  of  d  is 
increased  both  the  thickness-to-chord  ratio,  the  leading- 
edge  radius,  and  the  maximum  lift  coefficient  increase. 
These  airfoils  are  summarized  in  Table  5. 


TABLE  5 

T-SERIES  AIRFOIL  SUMMARY  (A  =  1/4) 


Airfoil 

d 

t/c  (%) 

T504 

0.5 

3.25 

0.0040 

0.605  /  5.5“ 

T604 

0.6 

4.3 

0.0067 

0.836  /  7.6“ 

T234 

0.666 

4.97 

0.0080 

0.943  /  8.6“ 

Case  5.  A  =  1/5.  For  the  values  of  d  evaluated  at  this 
value  of  A,  three  airfoils  meet  the  specified  objectives  for 
the  thickness-to-chord  ratio  and  the  leading-edge  radius. 

The  airfoil  profile,  leading-edge,  and  section  lift 
coefficient  for  Airfoils  T505,  T605  and  T6605  are  shown  in 
Figures  81,  82,  84;  85,  86,  88;  and  89,  90,  92, 
respectively.  These  figures  show  that  as  the  value  of  d  is 
increased  both  the  thickness-to-chord  ratio,  the  leading- 
edge  radius,  and  the  maximum  lift  coefficient  increase. 

These  airfoils  are  summarized  in  Table  6. 
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TABLE  6 


T-SERIES  AIRFOIL  SUMMARY  (A  =  1/5) 


Airfoil 

d 

t/c  (%) 

Cin«x  /  Stall  Angle 
(Re  =  3  X  10^) 

T505 

0.5 

2.89 

0.0032 

0.55  /  5.0° 

T605 

0.6 

3.82 

0.0055 

0.77  /  7.0° 

T6605 

0.666 

4.58 

0.0075 

0.903  /  8.3° 

Case  Comparisons 

Examining  the  pressure  distribution  and  Ci(a)  curves 
shows  the  significant  effect  the  leading-edge  radius  has  on 
the  characteristics  of  the  airfoil.  For  example,  focusing 
on  the  pressure  distributions  for  a  =  0®  in  Figures  33,  43, 
47,  and  51,  we  can  see  that  as  r^  increases  (which  increases 
as  these  figure  numbers  increase),  the  length  of  the  area 
with  a  favorable  pressure  gradient  near  the  leading-edge 
grows.  Also,  as  the  angle  of  attack  is  increased,  the 
leading-edge  radius  is  seen  to  have  a  great  effect  on  the 
pressure  distribution  near  the  leading-edge.  The  leading- 
edge's  influence  on  the  pressure  distribution  affects  the 
separation  of  the  flow,  because  of  the  adverse  pressure 
gradients.  The  results  of  this  separation  are  reflected  in 
Figures  40,  44,  48,  52,  56,  60,  64,  63,  72,  76,  80,  84,  88 
and  92.  The  maximum  lift  coefficient,  which  indicates  the 
onset  of  stall,  increases  as  the  leading-edge  radius 
increases. 

Also,  Figures  40,  44,  48,  52,  56,  60,  64,  68,  72,  76, 
80,  84,  88  and  92  show  the  maximum  lift  coefficient 
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increases  as  the  Reynolds  number  increases.  Reynolds  number 
is  basically  a  ratio  of  the  inertial  forces  to  the  viscous 
forces.  If  Reynolds  number  is  increased,  the  flow  is 
energized,  and  the  separation  length  along  the  airfoil  is 
reduced.  Also,  this  allows  for  a  higher  angle  of  attack 
before  ultimate  separation.  For  these  thin  airfoils,  the 
stall  characteristic  is  very  abrupt.  Therefore,  the  results 
we  see  are  most  significantly  influenced  by  the  latter  of 
the  two  effects. 

Design  Implications 

Using  the  power  law  relationship  was  very  successful 
for  developing  symmetric,  low-observable  general 
configuration  symmetric  airfoils  within  the  specified 
thickness-to-chord  and  leading-edge  radius  limits.  Figure 
93  shows  the  utilization  of  Equation  (67),  produces  input  a* 
distributions  resulting  in  airfoils  falling  into  the 
prescribed  t/c  band  for  0.25  £  d  <  0.75  for  various  values 
of  A.  In  the  presentation  of  Cases  1-5,  only  the  airfoils 
within  the  specified  objectives  for  thickness-to-chord  ratio 
and  leading-edge  radius  were  shown.  However,  Figure  94 
shows  the  results  for  all  the  values  of  d  evaluated  in  Cases 
1-5.  Figure  94  indicates  that  a  significant  portion  of  the 
airfoils  created  with  Equation  (67)  fall  simultaneously 
within  both  the  defined  t/c  and  r^  bands. 


The  most  significant  result  from  Figures  93  and  94  is 
that  the  designer  now  has  a  new  tool  for  developing  general 


low-observable  configuration,  symmetrical  airfoils.  Once 
the  designer  assesses  the  threat  that  an  airfoil  will  be 
facing,  a  determination  can  be  made,  based  on  the  threat 
frequency  as  discussed  earlier,  on  the  constraints  for  t/c 
and  rf  From  Figures  93  and  94,  the  designer  can  determine 
the  values  of  A  and  d  for  Equation  (67)  needed  to  get  the 
desired  t/c  and  r^.  However,  Figures  93  and  94  also  show 
the  designer  is  constrained,  for  a  given  r^,  to  a  t/c  band 
of  approximately  one  percent.  For  example,  if  the  designer 
determines  the  desired  airfoil  will  require  a  r^  =  0.0050, 
Figure  94  shows  that  an  airfoil  can  be  created  with  a 
thickness-to-cord  ratio  between  approximately  3.5  -  4.5 
percent  from  the  curves  shown. 

Comparison  of  T-Series  Airfoils  with  NACA  4-Diait  Airfoils 
The  T-Series  airfoils  need  to  be  evaluated  against  a 
similar  airfoil  geometry  to  contrast  their  performance 
characteristics.  One  of  the  major  design  parameters  for  the 
symmetric  T-series  airfoils  has  been  thickness-to-chord 
ratio.  The  need  for  symmetric  airfoils  based  on  thickness- 
to-chord  ratio  lends  itself  directly  to  symmetric  NACA  4- 
digit  airfoils  for  comparison.  The  NACA  4-digit  airfoils 
used  in  this  investigation  are  not  a  reference  to  common 
NACA  4-digit  airfoils  (i.e.  NACA  0006).  The  following 
symmetric  NACA  4-digit  airfoils  are  generated  from  the 
general  equation  for  the  thickness  distribution  for  NACA  4- 
digit  airfoils  given  by  Abbott  (10)  as 
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±y^.  =  Us  (  .2969y[x  -  .  126x  -  .3516x2  +  .2843x2  _  ^  loiSx^) 

The  program  contains  the  option  to  analyze  NACA  4-digit 
airfoils  with  the  potential-flow  analysis  method.  The 
procedure  requires  the  following  information  about  the 
airfoil;  the  amount  of  camber,  the  location  of  the  maximum 
camber,  and  the  thickness-to-chord  ratio.  This  information 
is  readily  available  since  these  three  pieces  of  information 
match  that  used  in  the  numbering  system  of  a  NACA  4-digit 
airfoil.  Eppler  (7)  addresses  the  more  specific  details  of 
this  option. 

An  airfoil  from  each  of  the  five  cases  are  chosen  for 
comparison  with  the  NACA  4-digit  airfoil.  Figures  95-101 
show  the  comparisons  of  the  airfoil  shape,  pressure 
distribution,  lift  coefficient,  lift  versus  drag,  lift-drag 
ratio  versus  angle  of  attack,  and  boundary-layer  development 
between  Airfoil  T3  and  the  corresponding  NACA  4-digit 
airfoil.  The  same  comparisons  are  shown  for  Airfoil  T22  in 
Figures  102-108,  for  Airfoil  T603  in  Figures  109-115,  for 
Airfoil  T234  in  Figures  116-122,  and  for  Airfoil  T605  in 
Figures  123-129. 

Figures  95,  102,  109,  116  and  123  show  the  T-series 
airfoils  have  a  thicker  front  and  more  tapered  tail  than  the 
NACA  4-digit  airfoils.  Also,  the  leading-edge  radius  is 
smaller  for  the  NACA  4-digit  airfoil.  Abbott  (10)  gives  the 
equation  for  the  leading-edge  radius  of  the  NACA  4-digit 
airfoils  to  be 
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(69) 


=  1.10l9(t/c)2. 


A  comparison  of  the  leading-edge  radii  is  summarized  in 
Table  7. 


TABLE  7 


LEADING-EDGE  RADII  (T-SERIES  &  NACA  4-DIGIT  AIRFOILS) 


T-Series 

Airfoil 

t/c 

(%) 

NACA  4-Digit 
Airfoil 

T3 

4.09 

0.0040 

0.0018 

T22 

4.75 

0.0060 

0.0025 

T603 

5.08 

0.0085 

0.0028 

T234 

4.97 

0.0080 

0.0027 

T605 

3.82 

0.0055 

0.0016 

Table  7  shows  that  the  leading-edge  radius  is  larger  for  the 
T-series  airfoils  than  the  NACA  4-digit  airfoil  with  the 
scime  t/c.  Therefore,  a  NACA  4-digit  airfoil  with  the  same 
leading-edge  radius  as  the  T-series  airfoil  would  have  a 
much  larger  thickness-to-chord  ratio.  These  thickness-to- 
chord  ratios  would  not  be  acceptable  to  the  specified 
constraints  of  this  investigation.  For  example,  in  Table  7, 
Airfoil  T3  has  the  smallest  value  of  r^.  Using  Equation 
(69),  a  NACA  4-digit  airfoil  with  the  same  value  of  r^  would 
result  in  t/c  =  6.03%. 

Figures  96,  103,  110,  117,  and  124  show  for  a  given 
angle  of  attack,  the  NACA  4-digit  airfoils  have  a  stronger 
adverse  pressure  gradient  on  the  upper  surface  than  the  T- 
Series  airfoils.  As  we  have  seen,  this  is  an  expected 
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result  for  an  airfoil  with  a  small  thickness-to-chord  ratio 


and  small  leading-edge  radius. 

The  T-Series  airfoils  have  larger  leading-edge  radii 
which  allows  for  a  smaller  adverse  pressure  gradient  in  the 
recovery.  Thus,  the  T-Series  airfoils  result  in  a  trade-off 
of  the  smallest  possible  leading-edge  radius  for  a  given 
adverse  pressure  gradient.  Figures  97,  104,  111,  118,  and 
125  show  a  benefit  of  this  tradeoff  in  the  section  lift 
coefficient.  The  benefit  of  this  tradeoff  for  the  five 
airfoils  are  summarized  in  Table  8. 


TABLE  8 

MAXIMUM  LIFT  COEFFICIENT  (T-SERIES  &  NACA  4-DIGIT  AIRFOILS) 


T-Series 

Airfoil 

t/C 

(%) 

^Itwx 

NACA  4-Digit 
Airfoil 

T3 

4.09 

0.77 

0.57 

T22 

4.75 

00 

CO 

• 

o 

0.68 

T603 

5.08 

0.95 

0.71 

T234 

4.97 

0.935 

0.71 

T605 

3.82 

0.77 

0.55 

Figures  97,  104,  111,  118,  and  125  and  Table  8  show  the  T- 
Series  airfoils  have  an  increased  stall  angle  of 
approximately  two  degrees  over  a  NACA  4-digit  airfoil  with 
the  same  t/c. 

An  increase  in  the  stall  angle  indicates  the  location 
of  the  separation  point  has  changed.  Figures  100,  101,  107, 
108,  114,  115,  121,  122,  126,  and  127  show  the  boundary- 
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layer  development  plots  for  the  upper  surfaces.  From 
Chapter  I,  the  significant  values  of  the  shape  function,  H32 
were  shown  to  be; 


H32  = 

1.51509 

(Laminar  separation) 

H32  = 

1.46 

(Turbulent  separation). 

H32  > 

1.58 

(Turbulent  reattachment). 

The  transition  line  on  these  figures  is  based  upon  Equation 
(50)  with  m  =  0  and  r  =  0.  It  must  be  mentioned  that,  in 
these  figures,  the  specific  point  where  laminar  separation 
or  transition  occur  according  to  the  analysis  may  not  be 
shown  intersecting  the  displayed  lines  for  these  events. 

This  is  a  problem  with  the  version  of  the  code.  In  this 
code,  the  integral  boundary-layer  method  evaluates  the 
criteria  between  the  points  plotted  in  the  figures. 

However,  the  results  of  the  evaluation  of  these  interior 
points  are  not  shown  in  the  figures. 

All  of  the  boundary-layer  development  analysis  was  done 
assuming  a  Re  =  3.0  x  10®.  Based  on  experiment,  this 
Reynolds  number  is  large  enough  so  Reynolds  number  effects 
are  small  (11).  In  Figures  100,  101,  107,  108,  114,  115, 

121,  122,  128,  and  129,  three  angles  of  attack  are  plotted 
on  each  of  the  boundary-layer  development  plots.  Each  plot 
shows  the  case  for  a  =  0.  The  other  two  angles  of  attack 
were  chosen  so  as  to  bracket  the  separation  point.  Figures 
100,  101,  107,  108,  114,  115,  121,  122,  128,  and  129  show 
that  the  boundary-layer  development  profile  for  a  =  0,  for 
both  the  T-Series  and  the  NACA  4-digit  airfoils,  indicate  a 


46 


transition  to  a  turbulent  boundary-layer. 

At  the  pre-separation  angle  of  attack  on  each  boundary- 
layer  development  plot  for  the  T-Series  airfoils,  Figures 
100,  107,  114,  121,  and  128  all  indicate  laminar  separation 
and  reattachment  with  a  turbulent  boundary-layer.  For  the 
NACA  4-digit  airfoils  corresponding  to  Airfoils  T3  and  T605, 
Figures  101  and  129  show  laminar  separation  followed  by  a 
reattachment  with  a  turbulent  boundary-layer.  However,  for 
the  NACA  4-digit  airfoils  corresponding  to  Airfoils  T22, 
T603,  and  T234,  Figures  108,  115,  and  122  show  transition 
before  laminar  separation  and  then  a  reattachment  of  a 
turbulent  boundary- layer. 

Figures  98,  105,  112,  119,  and  126  show  the  affect  of 
the  above  phenomenon  in  versus  plots.  The  transition 
to  a  turbulent  boundary-layer  before  laminar  separation 
reduces  the  drag  (4).  Figures  105,  112,  and  119  show  that 
when  this  occurs,  on  the  NACA  4-digit  airfoil,  the  bottom  of 
the  drag  bucket  is  lower  than  for  the  T-Series  airfoils. 

For  the  Airfoil  T3  and  T605  comparisons.  Figures  98  and  126 
show  that  there  is  no  significant  difference  in  the  bottom 
of  the  curve. 

Finally,  Figures  99,  106,  113,  120,  and  127  show  the 
comparison  of  the  lift-drag  ratio  curves.  Clancy  (12) 
indicates  this  relationship  is  important  because  it  can  be 
viewed  as  a  measure  of  airfoil  efficiency.  Overall,  there 
is  not  a  significant  difference  in  the  efficiency  of  T- 
Series  airfoils  versus  the  NACA  4-digit  airfoils. 
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V.  Conclusions  and  Recommendations 


Conclusions 

This  study  has  shown  the  applicability  of  PROFILE  in 
the  design  of  symmetric,  low-speed,  low-observable 
configuration  airfoils.  An  objective  of  this  investigation 
was  to  determine  any  limitations  of  the  program  for  this 
application.  As  discussed,  airfoils  with  small  leading-edge 
radii  and  small  thickness-to-chord  ratios  are  prone  to  early 
separation  and  abrupt  stall  characteristics.  For  airfoils 
with  very  small  leading-radii,  the  viscous  correction  for 
the  section  lift  coefficient  was  not  suitable  after  the 
stall  since  it  resulted  in  unrealistic  negative  values.  For 
the  T-Series  airfoils,  the  stall  angles  were  larger,  but  the 
same  effect  was  seen  in  the  inability  to  predict  post  stall 
lift  coefficients.  The  best  results  in  obtaining  this  class 
of  airfoils  was  achieved  if  the  secondary  pressure  recovery 
was  neglected  by  choosing  =  0  for  both  the  upper  and 
lower  surfaces. 

The  power  law  relationship  between  a*  and  x/c  based  on 
the  physical  relationships  between  the  pressure 
distribution,  a*,  leading-edge  shape,  and  the  thickness  of 
the  airfoil  proved  very  successful.  In  all  five  case 
studies  performed  for  the  different  coefficients  of  the 
power  law  relationship,  symmetrical  airfoils  were  produced 
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within  the  imposed  constraints  on  the  leading-edge  radius 
and  the  thickness-to-chord  radius.  The  curves  presented  in 
Figures  91  and  92  give  the  designer  a  new  ability  to 
determine  the  input  parameters  needed  to  produce  an  airfoil 
with  the  desired  thickness-to-chord  ratio  and  leading-edge 
radius.  However,  for  a  given  leading-edge  radius,  the 
designer  is  constrained  to  a  band  of  available  thickness-to- 
chord  ratios  with  a  width  of  approximately  one  percent. 

The  T-Series  airfoils,  produced  with  the  power  law 
relationship,  showed  improved  performance  compared  to  the 
NACA  4-digit  airfoils.  The  T-Series  airfoils  had  larger 
leading-edge  radii  than  an  equivalent  thickness  NACA  4-digit 
airfoils,  which  resulted  in  slightly  higher  drag.  However, 
the  stall  angles  were  increased  approximately  two  degrees. 
Therefore,  the  T-Series  airfoils  showed  improved  stall 
characteristics,  for  an  airfoil  meeting  the  low- 
observability  criteria,  compared  to  equivalent  thickness 
NACA  4-digit  airfoil . 

Recommendations  for  Further  Study 

1.  Investigate  the  possibility  of  extending  the  power 
law  relationship  for  thin  cambered  airfoils. 

2.  Conduct  a  parcunetric  study  in  which  input 

parameters  held  constant  in  this  investigation,  are  varied. 
The  closure  contribution  parameters  (pg,  Kj^,  and  could 

be  varied  to  determine  their  effects  on  the  airfoil 
thickness  and  aerodynamic  characteristics  due  to  the  change 
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in  the  shape  of  the  t railing-edge. 

3.  Investigate  the  performance  of  airfoils  created 
with  a  transonic  analysis  code. 
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Figure  2.  Arc  Limit  Description  (1:129) 
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Figure  3.  Pressure  Recovery  and  Closure  Contribution  (1:128) 
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and  Pressure  Distributions  -  t/c  =  14.27%  -  Incompressible 


igure  5.  Sensitivity  of  Pressure  Distribution  to  Changes  in  Location 
of  Specified  Main  Pressure  Recovery  -  Airfoil  G48  -  a  =  2° 
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Coefficient  to  the  Location  of  the  Secondary  Pressure  Recovery  -  Airfoil  G48 
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Figure  8.  Velocity  Distributions  -  Airfoil  SH9  -  Incompressible  - 


Figure  9.  Velocity  Distributions  -  Airfoil  SH9  -  Incompressible  - 


Without  B.  L.  displacement  iteration 
With  B.  L.  displacement  iteration 


Figure  10.  Pressure  Distributions  -  NACA  4415  - 


Figure 


14.  Airfoil  G31  Leading-Edge  -  t/c  =  4.95% 


Figure  15.  Airfoil  G31  -  C,  vs  a  Curve 


ickness-to-Chord  Ratio 


Figure  18.  Pressure  Distributions  -  NACA  0006  vs  NACA  0012  -  Incompressible 


Figure  19.  TA2-TA5  Curves  -  Proposed  a  -X  Relationships 


Figure  20.  Airfoil  TA2  -  t/c  =  12.5% 


21.  Airfoil  TA2  Leading-Edge 


Figure  22.  Pressure  Distributions  -  TA2  -  Incompressible 


igure 


24.  Airfoil  TA3 
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Figure  25.  Pressure  Distributions  -  TA3  -  Incompressible 


Figure  28.  Pressure  Distributions  -  TA4  •  Incompressible 


5  Leading-Edge  -  t/c  =  2.37% 
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Figure  31.  Pressure  Distributions  -  TA5  -  Incompressible 
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Figure  32.  Power  Law  Relationships  {A=1) 
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Figure  33.  Power  La 
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Figure  39.  Pressure  Distributions  -  T25  -  Incompressible 
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Figure  40.  Airfoil  T25  -  C,,  vs  a  Curves 


ading-Edge  -  t/c  =  3.59% 
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Figure  43.  Pressure  Distributions  -  T30  -  Incompressible 
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Airfoil  T3  Leading-Edge  -  t/c  =  4.09% 
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Figure  47.  Pressure  Distributions  -  T3  -  Incompressible 
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Section  angle  of  attack,  a  (deg.) 
Figure  48.  Airfoil  T3  -  C,,  vs  a  Curves 


igure 


Leading-Edge  -  t/c  =  5.17% 
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Figure  51.  Pressure  Distributions  -  T4  -  Incompressible 
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T402  -  t/c  =  3.44% 


54.  Airfoil  T402  Leading-Edge  -  t/c  =  3.44% 


Figure  55.  Pressure  Distributions  -  T402  -  Incompressible 
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Figure  56.  Airfoil  T402  -  C|,C^  vs  a  Curves 
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igure  57.  Airfoil  T22  -  t/c  =  4.75% 


e  58.  Airfoil  T22  Leading-Edge  -  t/c  =  4.75% 


Figure  59.  Pressure  Distributions  -  T22  -  Incompressible 
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Figure  60.  Airfoil  T22  -  C,.  vs  a  Curves 


62.  Airfoil  T503  Leading-Edge  -  t/c  =  3.79% 
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Figure  63.  Pressure  Distributions  -  T503  -  Incompressible 
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Figure  64.  Airfoil  T503  -  C,,  vs  a  Curves 


T603  -  t/c  =  5.08% 
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Figure  67.  Pressure  Distributions  -  T603  -  Incompressible 
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Section  angle  of  attack,  a  (deg.) 
Figure  68.  Airfoil  T603  -  C,,  vs  a  Curves 


Figure  69.  Ai 
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Figure  71.  Pressure  Distributions  -  T504  -  Incompressible 
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Figure  72.  Airfoil  T504  -  C,,  vs  a  Curves 


Figure  73.  Airfoil  T604  -  t/c  =  4.3% 
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Figure  75.  Pressure  Distributions  -  T604  -  Incompressible 
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Section  angle  of  attack,  a  (deg.) 

Figure  76.  Airfoil  T604  -  C,,  vs  a  Curves 


igure 


Figure  79.  Pressure  Distributions  -  T234  -  Incompressible 
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Figure  81.  Airfoil  T505  -  t/c  =  2.89 
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Figure  83.  Pressure  Distributions  -  T505  -  Incompressible 
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Figure  84.  Airfoil  T505  -  C,.  vs  a  Curves 


0.010 


138 


139 


Figure  87.  Pressure  Distributions  -  T605  -  Incompressible 
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Section  angle  of  attack,  a  (deg.) 
Figure  88.  Airfoil  T605  -  C,,  vs  a  Curves 
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Figure  89.  Airfoil  T6605  -  t/c  =  4.58% 


Leading-Edge  -  t/c  =  4.58% 
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Figure  91 .  Pressure  Distributions  -  T6605  -  Incompressible 
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Figure  92.  Airfoil  T6605  -  C,,  vs  a  Curves 
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Figure  94.  Leading-Edge  Radius  vs  Thickness-to-Chord  Ratio 
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Figure  103.  Airfoil  T22  vs  NACA  4-Digit  Airfoil  Pressure  Distribution 

a  =  6°  -  Incompressible 
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Figure  104.  C,(a)  Curve  -  Airfoil  T22  vs  NACA  4-Digit  Airfoil 
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Figure  105.  C,  -  plot  -  Airfoil  T22  vs  NACA  4-cliglt  airfoil 

t/c  =  4.75% 


Re  =  3.0  X  10' 
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Figure  110.  Airfoil  T603  vs  NACA  4-Digit  Airfoil  Pressure  Distribution 

a  =  6°  -  Incompressible 
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-16  -8  0  8  16 
Section  angle  of  attack,  a  (deg.) 

Figure  111.  C|(a)  Curve  -  Airfoil  T603  vs  NACA  4-Digit  Airfoil 
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Figure  117.  Airfoil  T234  vs  NACA  4-Digit  Airfoil  Pressure  Distribution 

a  =  6®  -  Incompressible 
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Figure  ViB.  C|(a)  Curve  -  Airfoil  T234  vs  NACA  4-Digit  Airfoil 
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Re  =  3.0  X  10' 


Figure  125.  C,(a)  Curve  -  Airfoil  T605  vs  NACA  4-Digit  Airfoil 
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Figure  126.  C,  -  plot  -  Airfoil  T605  vs  NACA  4-cligit  airfoil 

t/c  =  3.82% 
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Appendix  A;  Summary  of  the  Pairs 

This  appendix  contains  the  pairs  determined 

for  the  TA-series  and  the  T-series  airfoils.  The  TA-series 
airfoils  were  obtained  by  the  hand-drawn  curves  of 
decreasing  a*.  The  T-series  airfoils  were  obtained  for  the 
power  law  relationship  defined  by  Equation  (67). 
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TABLE  A.l 


SUMMARY  OF  (v^,a^*)  PAIRS  FOR  TA-SERIES  AIRFOILS 
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TABLE  A. 2 


SUMMARY  OF  PAIRS  FOR  T-SERIES  AIRFOILS  (A  =  1) 


5 


10 


15 


20 


25 


26 


27 


28 


28.3 


28.7 


29 


29.3 


29.6 


1.02 

1.02 

1.07 

1.09 

1.18 

1.23 

1.41 

1.51 

1.96 

2.25 

2.19 

2.56 

2.53 

3.04 

3.05 

3.81 

3.29 

4.18 

3.67 

4.76 

4.41 

5.45 

4.47 

6.03 

00 

o 

• 

in 

7.03 

6.68 

9.78 

-6.68 

-9.78 

-5.08 

-7.03 

-4.47 

-6.03 

-4.41 

-5.45 

-3.67 

-4.76 

-3.29 

-4.18 

-3.05 

-3.81 

-2.53 

-3.04 

-2.19 

-2.56 

-1.96 

-2.25 

-1.41 

-1.51 

-1.18 

-1.23 

-1.07 

-1.09 

-1.02 

-1.02 

1.02 


1.26 


1.58 


2.45 


2.48 


3.44 


4.55 


4.89 


5.65 


6.57 


7.35 


8.7 


12.57 


-12.57 


-8. 7 


-7.35 


-6.57 


-5.65 


-4.89 


-4.55 


-3.44 


-2.84 


-2.45 


-1.58 


-1.26 


-1.02 


1.028 


1.12 


1.32 


1.74 


2.93 


3.5 


8.0 


9.6 


10.9 


13.4 


20.9 


-20.9 


-13.4 


-10.9 


-9.6 


-8.0 


-3.5 


-2.93 


-1.74 


-1.32 


-1.12 


-1.028 
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A. 3 


l85 


TABLE  A.  4 


SUMMARY  OF  PAIRS  FOR  T-SERIES  AIRFOILS  (A  =  1/3 


i 

''i 

T503 

* 

T603 

* 

1 

5 

.34 

.34 

2 

10 

.38 

.39 

3 

15 

.47 

.  .5 

4 

20 

.66 

.76 

5 

25 

1.3 

1.69 

6 

26 

1.6 

2.2 

7 

27 

2.1 

3.1 

6 

28 

3.2 

5.1 

9 

28.3 

3.7 

5.8 

10 

28.7 

4.5 

7.5 

11 

29 

5.6 

9.9 

12 

29.3 

6.67 

12.1 

13 

29.6 

8.6 

16.5 

14 

0 

14.9 

31.9 

15 

30.4 

-14.9 

-31.9 

16 

30.7 

-8.6 

-16.5 

17 

31 

-6.67 

-12.1 

18 

31.3 

-5.6 

-9.9 

19 

31.7 

-4.5 

-7.5 

20 

32 

-3.7 

-5.8 

21 

33 

-3.2 

-5.1 

22 

34 

-2.1 

-3.1 

23 

35 

-1.6 

-2.2 

24 

-1.3 

-1.69 

25 

-.66 

-.76 

26 

50 

-.47 

-.5 

55 

-.38 

-.39 

60 

-.34 

-.34 
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A. 6 


L88 


The  purpose  of  this  appendix  is  to  explain  the  function 
of  the  major  input  lines  used  in  this  investigation.  Eppler 
(7)  details  all  the  input  options  for  this  program. 

The  sequence  of  the  program's  execution  is  managed  by 
the  named  input  lines.  Once  the  program  reads  the  name,  it 
goes  to  the  correct  section  of  the  code.  After  completing 
the  assigned  task,  it  returns  back  to  the  main  program  and 
reads  the  next  input  line.  The  following  is  a  sample  input 
file: 

REMOl 

TRA21  1000  -1  -.7  -.5  -.1  0  .1  .5  .7  1  1000 

TRAl  22.5  0  0  2  37.5  -2  60  0 

TRA21  4  29  1  .1  .95  4  29  1  .1  . 95  3  .302 

PAN  13 

ALFA213  7  -6 . , -4 . , -2 . , 0 . , 2 . , 4 . , 6 . 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 

ENDE 

The  REMOl-line  is  used  here  to  toggle  between  the  free 
and  formatted  modes  of  input.  The  character  ”1"  in  the  5th 
position  toggles  to  the  free  format  mode. 

The  TRAl-  and  TRA2-line  contain  the  information  needed 
in  the  airfoil  design.  The  first  TRA2 1-line  specifies  the 
information  for  putting  additional  points  at  the  leading 
edge.  This  information  is  triggered  into  use  by  specifying 
the  character  ”2"  as  the  14th  character  in  the  second  TRA21 
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line.  The  TRAl-line  used  here  contains  the  (Vi,a^*)  pairs, 
where  is  the  specification  of  the  i-th  arc  segment 
related  to  by 


V 


i 


N 


360 


(70) 


The  second  TRA21-line  specifies  the  following  (rotating 
input  line  into  vertical ) : 

TRA21 

4  A.*  =  beginning  of  closure  contribution,  upper  side 

29  A.  =  beginning  of  main-pressure  region,  upper  side 

1  RMSyg  =  recovery  specification  mode,  upper  side 
(determines  interpretation  of  next  two 


inputs ) 

RMSyg  =  0 

1 

2 

3 

.1  =  K 

w' 

y- 

y 

.95  =  p 

w 

o 

o' 

4  I*  =  beginning  of  closure  contribution,  lower  side 

29  J.  =  beginning  of  main-pressure  region,  lower  side 

1  RMSj^g  =  recovery  specification  mode,  lower  side 
(determines  interpretation  of  next  two 
inputs ) 

RMSia  =0  12  3 

.1  =  K  S'  fi  Ji 

.95  =  p  S  5  5' 

3  ITMOD  =  trailing  edge  iteration  mode 
.3  Kr  =  desired  value  of  Kg,  Ks  “ 

0  =  tolerance  on 
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2  Triggers  use  of  first  TRA2 1-line  for  additional 
points  at  the  leading  edge 

The  Pan-line  switches  the  program  from  the  design  mode 
to  the  analysis  mode  (panel  method).  The  ALFA-line 
specifies  the  a’s,  the  angles  of  attack  of  the  incoming  flow 
in  degrees,  at  which  the  airfoil  should  be  evaluated.  The 
RE-line  specifies  the  Reynolds  numbers  for  which  the  airfoil 
should  be  evaluated  and  which  transition  mode  to  use  in  the 
boundary-layer  analysis.  The  ENDE-line  terminates  the 
program's  execution. 
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Appendix  C:  Summary  of  Input  Files 
This  appendix  contains  the  input  files  used  to  create 
the  airfoils  presented  in  this  thesis. 
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Airfoil  TA2 


REN01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 
TRA1  5  0  10  2.2  15  5.5  20  9  25  13  26  19  27  19.6  28  21.1 

TRA1  28.3  23  28.7  24  29  25.1  29.3  26  29.6  27.3 

TRA1  0  30  30.4  -30  30.7  -27.3  31  -26  31.3  -25.1 

TRA1  31.7  -24  32  -23  33  -21.1  34  -19.6  35  -19  40  -13 

TRA1  45  -9  50  -5.5  55  -2.2  60  0 

TRA21  401  .1  .9401  .1  .9  3  .3  0  2 
PAN  13 

ALFA213  13  -12., -10., -8.. -6., -4., -2., 0., 2.. 4. ,6. ,8. ,10. ,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENOE 


Airfpil  TA3 
REH01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 
TRA1  5  0  10  .9  15  1.7  20  3  25  5  26  9.5  27  11  28  13 

TRA1  28.3  16  28.7  16.8  29  18.2  29.3  19.7  29.6  21.9 

TRA1  0  30  30.4  -30  30.7  -21.9  31  -19.7  31.3  -18.2 

TRA1  31.7  -16.8  32  -16  33  -13  34  -11  35  -9.5  40  -5 

TRA1  45  -3  50  -1.7  55  -.9  60  0 

TRA21  4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12.,-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 


Airfoil  TA4 
REN01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRA1  5  0  10  .2  15  .7  20  .9  25  1  26  2.5  27  3.8  28  5.5 

TRA1  28.3  9.7  28.7  10.8  29  13.2  29.3  16  29.6  20 

TRA1  0  25  30.4  -25  30.7  -20  31  -16  31.3  -13.2 

TRA1  31.7  -10.8  32  -9.7  33  -5.5  34  -3.8  35  -2.5  40  -1 

TRA1  45  -.9  50  -.7  55  -.2  60  0 

TRA21  401  .1  .9401  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12.,-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 

Airfoil  TA5 
REN01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRA1  5  0  10  0  15  0  20  0  25  0  26  0  27  .01  28  .2 

TRA1  28.3  1.2  28.7  1.9  29  3  29.3  4.5  29.6  7.8 

TRA1  0  30  30.4  -30  30.7  -7.8  31  -4.5  31.3  -3 

TRA1  31.7  -1.9  32  -1.2  33  -.2  34  -.01  35  0  40  0 

TRA1  45  0  50  0  55  0  60  0 

TRA21  401  .1  .95  4  0  1  .1  .95  3  .3  0  2 

PAN  13 

ALFA213  13  -12.,-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENOE 
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Airfoil  125 


REH01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRA1  5  1.02  10  1.07  15  1.18  20  1.41  25  1.96  26  2.19  27  2.53  28  3.05 

TRA1  28.3  3.29  28.7  3.67  29  4.41  29.3  4.47  29.6  5.08 

TRA1  0  6.68  30.4  -6.68  30.7  -5.08  31  -4.47  31.3  -4.41 

TRA1  31.7  -3.67  32  -3.29  33  -3.05  34  -2.53  35  -2.19  40  -1.96 

TRA1  45  -1.41  50  -1.18  55  -1.07  60  -1.02 

TRA21  4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12. ,-10. , -8. , -6. , -4. , -2. .0. ,2. ,4. .6. ,8. .10. ,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 


Airfoil  130 

REN01 
TRA21 
TRAI 
TRA1 
TRAI 
TRAI 
TRAI 
TRA21 
PAN  13 
ALFA213  13 
RE  130  0 
ENDE 

Airfoil  T3 

REH01 
TRA21 
TRAI 
TRA1 
TRAI 
TRAI 
TRAI 
TRA21 
PAN  13 
ALFA213  13 
RE  130  0 
ENDE 

Airfoil  T4 
REN01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRAI  5  1.028  10  1.12  15  1.32  20  1.74  25  2.93  26  3.5  27  4.4 

TRAI  28  6  28.3  6.7  28.7  8  29  9.6  29.3  10.9  29.6  13.4 

TRAI  0  20.9  30.4  -20.9  30.7  -13.4  31  -10.9  31.3  -9.6 

TRAI  31.7  -8  32  -6.7  33  -6  34  -4.4  35  -3.5  40  -2.93 

TRAI  45  -1.74  50  -1.32  55  -1.12  60  -1.028 

TRA21  4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12. ,-10. , -8. , -6. , -4. , -2. ,0. ,2. ,4. ,6. ,8. ,10. ,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 


1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 
5  1.02  10  1.1  15  1.26  20  1.58  25  2.45  26  2.84  27  3.44 
28  4.55  28.3  4.89  28.7  5.65  29  6.57  29.3  7.35  29.6  8.7 
0  12.57  30.4  -12.57  30.7  -8.7  31  -7.35  31.3  -6.57 
31.7  -5.65  32  -4.89  33  -4.55  34  -3.44  35  -2.84  40  -2.45 
45  -1.58  50  -1.26  55  -1.1  60  -1.02 
4  0  1  .1  .9  4  0  1  .1  .9  3  0  0  2 

-12. ,-10. ,-8. ,-6. ,-4. ,-2. ,0. ,2., 4., 6., 8. ,10. ,12. 

03.0  3000  3.0  6000  3.0  9000 


1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

5  1.02  10  1.09  15  1.23  20  1.51  25  2.25  26  2.56  27  3.04  28  3.81 

28.3  4.18  28.7  4.76  29  5.45  29.3  6.03  29.6  7.03 

0  9.78  30.4  -9.78  30.7  -7.03  31  -6.03  31.3  -5.45 

31.7  -4.76  32  -4.18  33  -3.81  34  -3.04  35  -2.56  40  -2.25 

45  -1.51  50  -1.23  55  -1.09  60  -1.02 

4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

-12., -10., -8., -6., -4., -2., 0., 2., 4., 6. ,8., 10., 12. 

03.0  3000  3.0  6000  3.0  9000 
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Airfoil  TAOe 


REH01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 
TRAI  5  .51  10  .56  15  .66  20  .86  25  1.5  26  1.7  27  2.2  28  3 

TRA1  28.3  3.4  28.7  4  29  4.8  29.3  5.4  29.6  6.7 

TRAI  0  10.4  30.4  -10.4  30.7  -6.7  31  -5.4  31.3  -4.8 

TRAI  31.7  -4  32  -3.4  33  -3  34  -2.2  35  -1.7  40  -1.5 

TRAI  45  -.86  50  -.66  55  -.56  60  -.51 

TRA21  4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12., -10. , -8., -6., -4., -2., 0., 2.. 4. .6. ,8.. 10., 12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 

Airfoil  T22 

REHOI 
TRA21 
TRAI 
TRAI 
TRAI 
TRAI 
TRAI 
TRA21 
PAN  13 
ALFA213 
RE  130 
ENOE 

Airfoil  TSOS 

REHOI 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRAI  5  .34  10  .38  15  .47  20  .66  25  1.3  26  1.6  27  2.1  28  3.2 

TRAI  28.3  3.7  28.7  4.5  29  5.6  29.3  6.67  29.6  8.6 

TRAI  0  14.9  30.4  -14.9  30.7  -8.6  31  -6.67  31.3  -5.6 

TRAI  31.7  -4.5  32  -3.7  33  -3.2  34  -2.1  35  -1.6  40  -1.3 

TRAI  45  -.66  50  -.47  55  -.38  60  -.34 

TRA21  4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12.,-10.,-8..-6.,-4.,-2.,0.,2.,4.,6.,8.,10.,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 

Airfoil  T6CS 
REN01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRAI  5  .34  10  .39  15  .5  20  .76  25  1.69  26  2.2  27  3.1  28  5.1 

TRAI  28.3  5.8  28.7  7.5  29  9.9  29.3  12.1  29.6  16.5 

TRAI  0  31.9  30.4  -31.9  30.7  -16.5  31  -12.1  31.3  -9.9 

TRAI  31.7  -7.5  32  -5.8  33  -5.1  34  -3.1  35  -2.2  40  -1.69 

TRAI  45  -.76  50  -.5  55  -.39  60  -.34 

TRA21  401  .1  .9401  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12.,-10.,-8..-6.,-4.,-2.,0.,2.,4.,6.,8.,10.,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 


1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 
5  .51  10  .58  15  .7  20  1  25  1.9  26  2.4  27  3.2  28  4.65 
28.3  5.4  28.7  6.75  29  8.45  29.3  10  29.6  12.9 
0  22.35  30.4  -22.35  30.7  -12.9  31  -10  31.3  -8.45 
31.7  -6.75  32  -5.4  33  -4.65  34  -3.2  35  -2.4  40  -1.9 
45  -1  50  -.7  55  -.58  60  -.51 
4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

13  -12.. -10., -8., -6., -4., -2., 0., 2., 4., 6.. 8.. 10., 12. 

0  03.0  3000  3.0  6000  3.0  9000 
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Airfoil  T5M 
REH01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

THAI  5  .26  10  .29  15  .35  20  .5  25  .96  26  1.2  27  1.6  28  2.4 

TRAI  28.3  2.7  28.7  3.4  29  4.2  29.3  5  29.6  6.4 

TRA1  0  11.2  30.4  -11.2  30.7  -6.4  31  -5  31.3  -4.2 

TRAI  31.7  -3.4  32  -2.7  33  -2.4  34  -1.6  35  -1.2  40  -.96 

TRAI  45  -.5  50  -.35  55  -.29  60  -.26 

TRA21  4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12., -10. ,-8.. -6., -4., -2., 0., 2., 4.. 6. .8. ,10.. 12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 


Airfoil  T604 
REH01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRAI  5  .26  10  .29  15  .38  20  .57  25  1.27  26  1.65  27  2.3  28  3.6 

TRAI  28.3  4.3  28.7  5.6  29  7.4  29.3  9.1  29.6  12.4 

TRAI  0  23.9  30.4  -23.9  30.7  -12.4  31  -9.1  31.3  -7.4 

TRAI  31.7  -5.6  32  -4.3  33  -3.6  34  -2.3  35  -1.65  40  -1.27 

TRAI  45  -.57  50  -.38  55  -.29  60  -.26 

TRA21  401  .1  .9401  .1  .9  3  .3  0  2 
PAN  13 

ALFA213  13  -12.,-10.,-8.,-6.,-4.,-2.,0.,2.,4.,6.,8.,10.,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 


Airfoil  T234 
REH01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRAI  5  .26  10  .3  15  .39  20  1.12  25  1.5  26  2  27  2.95 

TRAI  28  4.9  28.3  6  28.7  8  29  10.8  29.3  13.5  29.6  19 

TRAI  0  39.48  30.4  -39.48  30.7  -19  31  -13.5  31.3  -10.8 

TRAI  31.7  -8  32  -6  33  -4.9  34  -2.95  35  -2  40  -1.5 

TRAI  45  -1.12  50  -.39  55  -.3  60  -.26 

TRA21  401  .1  .9401  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12. ,-10. , -8. , -6. , -4. , -2. ,0. ,2. ,4. ,6. ,8. ,10. ,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 

Airfoil  T505 
REH01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 
TRAI  5  .21  10  .23  15  .28  20  .4  25  .77  26  .97  27  1.29  28  1.9 

TRAI  28.3  2.2  28.7  2.8  29  3.38  29.3  4  29.6  5.16 

TRAI  0  8.94  30.4  -8.94  30.7  -5.16  31  -4  31.3  -3.38 

TRAI  31.7  -2.8  32  -2.2  33  -1.9  34  -1.29  35  -.97  40  -.77 

TRAI  45  -.4  50  -.28  55  -.23  60  -.21 

TRA21  401  .1  .9401  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12. ,-10. , -8. , -6. , -4. , -2. ,0. ,2. ,4. ,6. ,8., 10. ,12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 
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Airfoil  TAOS 
REM01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

THAI  5  .20  10  .23  15  .3  20  .46  25  1.  26  1.3  27  1.87  28  3. 

TRA1  28.3  3.5  28.7  4.5  29  5.6  29.3  7.2  29.6  9.9 

TRAI  0  19.1  30.4  -19.1  30.7  -9.9  31  -7.2  31.3  -5.6 

TRA1  31.7  -4.5  32  -3.5  33  -3.  34  -1.87  35  -1.3  40  -1. 

TRAI  45  -.46  50  -.3  55  -.23  60  -.20 

TRA21  401  .1  .9401  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12.,-10..-8.,-6.,-4.,-2.,0..2.,4..6..8.,10..12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENOE 


Airfoil  T6605 
REH01 

TRA21  1000  -1  -.7  -.5  -.3  -.2  -.1  0  .1  .2  .3  .5  .7  1  1000 

TRAI  5  .21  10  .24  15  .31  20  .5  25  1.2  26  1.62  27  2.4  28  3.9 

TRAI  28.3  4.8  28.7  6.4  29  8.6  29.3  10.8  29.6  15.2 

TRAI  0  31.6  30.4  -31.6  30.7  -15.2  31  -10.8  31.3  -8.6 

TRAI  31.7  -6.4  32  -4.8  33  -3.9  34  -2.4  35  -1.62  40  -1.2 

TRAI  45  -.5  50  -.31  55  -.24  60  -.21 

TRA21  4  0  1  .1  .9  4  0  1  .1  .9  3  .3  0  2 

PAN  13 

ALFA213  13  -12. , -10. , -8. , -6. , -4., -2. ,0., 2. ,4. ,6. ,8., 10.. 12. 

RE  130  0  03.0  3000  3.0  6000  3.0  9000 
ENDE 
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